改进的均衡滤波器在位场数据边界识别中的应用   总被引：9，自引：8，他引：1       下载免费PDF全文

马国庆  黄大年  于平  李丽丽 《地球物理学报》2012,55(12):4288-4295

边界识别是进行位场数据解释时一项必不可少的任务.现有的边界识别滤波器存在识别边界发散和深部地质体边界模糊的缺点.本文提出增强型均衡滤波器,可有效地改善上述缺点.该滤波器是利用不同阶导数之间的组合来进行地质体边界的识别,并在运算中引入一种计算高阶垂直导数的稳定算法.通过理论模型试验证明增强型均衡滤波器能使浅部与深部地质体的边界同时清晰地显示,且相对于其它边界识别滤波器能更加准确和清晰地识别出地质体的边界.最后将增强型均衡滤波器应用于实测位场数据的解释,根据其识别结果可容易地划分出断裂的水平位置及不同地层之间的界线,并能发现更多的细节信息.  相似文献   

重磁位场转换计算中迭代法的综合分析与研究   总被引：9，自引：6，他引：3       下载免费PDF全文

姚长利  李宏伟  郑元满  孟小红  张聿文 《地球物理学报》2012,55(6):2062-2078

处理转换计算在重磁资料解释中发挥着重要的作用,但一些计算如向下延拓、化极等有时是很不稳定的,在频率域中则表现为其转换因子具有明显的放大作用,所以其FFT理论计算结果是不稳定的.因此,很多研究工作都是围绕增加计算的稳定性、提高计算效果进行的,其中迭代法是近来在研究中受到普遍重视的方法技术,并取得了较好的成果.但也存在对迭代法研究还不够深入,对其存在的缺点认识不够充分、客观等问题,例如,迭代法进行延拓及化极等计算时,对一些具体应用虽能在一定程度上获得较好的计算结果,但却存在计算结果并不会随着迭代次数的增加而得到持续改善的问题,对于原本不稳定的计算,迭代法在迭代次数比较大时,所得的结果依然是不稳定的.为此,本文在对迭代法进行分析研究的基础上,进一步推导了迭代法的通式,并分析了对迭代法收敛性影响的各种因素.分析结果表明:迭代法收敛到FFT理论直接计算结果的决定因素是计算过程中如何选择原始数据到目标数据的映射函数;在选择了合适的映射函数的情况下,迭代次数不仅仅是决定计算成本,而是决定结果好坏的关键因素;增加迭代次数虽然能够使计算收敛到FFT直接计算理论结果,但如果该理论结果本身就是不稳定的,则迭代法计算如果收敛,也是收敛到一个不稳定的结果.所以针对位场处理转换中一些不稳定计算采用迭代法,并没有从根本上解决计算的不稳定性问题.  相似文献   

解释位场全张量数据的张量局部波数法及其与常规局部波数法的比较   总被引：6，自引：5，他引：1       下载免费PDF全文

马国庆  杜晓娟  李丽丽 《地球物理学报》2012,55(7):2450-2461

全张量探测技术以其信息量大、精度高、干扰小等优点在地球物理领域中得到广泛应用.本文提出采用张量局部波数法来进行位场全张量数据的解释,首先给出了张量局部波数的定义,然后推导出利用张量局部波数法进行反演的基本公式.本文方法在进行张量数据反演时无需事先知道场源体的类型(构造指数)即可获得场源体的位置信息,且可根据位置参数对场源体的类型进行估计.通过理论模型证明张量局部波数法可以很好地完成位场全张量数据的反演工作,并将其与常规局部波数法进行对比,证明全张量局部波数法的反演结果更加准确,即使在测点分布不合理的情况下,张量局部波数法仍可以获得准确的结果.最后应用张量局部波数法对美国得克萨斯州实测重力数据进行了反演,其反演结果与已有的研究成果相一致.  相似文献   

MCC转为带状MCSs过程中水平涡度的变化与暴雨的关系   总被引：4，自引：0，他引：4  

丁治英  高松  常越 《热带气象学报》2013,29(4):540-550

利用实况资料和WRF中尺度数值模式对2010年6月18—19日的一次MCC转带状MCSs的暴雨过程进行数值模拟与诊断分析。结果表明:850 hPa西南涡和切变线的形成与维持是影响此次暴雨产生的中尺度系统,前期MCC的形成到成熟以低涡降水为主,后期的圆形MCC转为带状MCSs主要为切变线降水。在雨区附近,u、v的垂直切变所形成的强水平涡度造成的旋转,对应垂直环流的上升支可触发暴雨产生,垂直方向上u、v不同的分布可形成不同的垂直环流。低涡与切变线附近的水平涡度有明显差异,这种差异导致暴雨形成的原因不同,低涡暴雨主要由v的垂直切变造成,切变线暴雨主要由u、v的垂直切变共同作用,本次过程中v的垂直切变构成了沿切变线的东西向雨带,u的垂直切变沿纬向的不均匀性引起的垂直运动与切变线上MCSs的生成、发展和多雨团的形成关系密切。低涡、切变线降水中心附近的正倾侧项（水平涡度向垂直正涡度转换）也有类似的差异,低涡的转换主要由?v/?p<0决定,切变线的转换主要由-?u/?p>0决定。水平涡度向垂直涡度的转换尺度较小,易在平均状态下被忽略。倾侧项主要有利于暴雨的加强,但对西南涡、切变线的发展贡献较小。   相似文献   

Multiridge Euler deconvolution   总被引：1，自引：0，他引：1  

Giovanni Florio  Maurizio Fedi 《Geophysical Prospecting》2014,62(2):333-351

Potential field interpretation can be carried out using multiscale methods. This class of methods analyses a multiscale data set, which is built by upward continuation of the original data to a number of altitudes conveniently chosen. Euler deconvolution can be cast into this multiscale environment by analysing data along ridges of potential fields, e.g., at those points along lines across scales where the field or its horizontal or vertical derivative respectively is zero. Previous work has shown that Euler equations are notably simplified along any of these ridges. Since a given anomaly may generate one or more ridges we describe in this paper how Euler deconvolution may be used to jointly invert data along all of them, so performing a multiridge Euler deconvolution. The method enjoys the stable and high‐resolution properties of multiscale methods, due to the composite upward continuation/vertical differentiation filter used. Such a physically‐based field transformation can have a positive effect on reducing both high‐wavenumber noise and interference or regional field effects. Multiridge Euler deconvolution can also be applied to the modulus of an analytic signal, gravity/magnetic gradient tensor components or Hilbert transform components. The advantages of using multiridge Euler deconvolution compared to single ridge Euler deconvolution include improved solution clustering, increased number of solutions, improvement of accuracy of the results obtainable from some types of ridges and greater ease in the selection of ridges to invert. The multiscale approach is particularly well suited to deal with non‐ideal sources. In these cases, our strategy is to find the optimal combination of upward continuation altitude range and data differentiation order, such that the field could be sensed as approximately homogeneous and then characterized by a structural index close to an integer value. This allows us to estimate depths related to the top or the centre of the structure.  相似文献   

基于时-空ETAS模型给出的川滇地区背景地震活动和强震潜在危险区   总被引：5，自引：2，他引：3       下载免费PDF全文

蒋长胜  庄建仓 《地球物理学报》2010,53(2):305-317

利用基于时-空传染型余震序列（Epidemic Type Aftershock Sequence, 简称ETAS）模型的随机除丛法,重新审视了2008年5月12日汶川M_S8.0地震前可能存在的长期地震活动异常,研究了川滇地区背景地震活动特征,并评估了当前的强震危险状态.对川滇地区1970年以来的M_L3.0以上的背景地震和丛集地震活动的研究结果表明,该地区地震丛集特征明显、时空分布很不均匀、地震序列常有前震事件.直接将概率值作为地震计数的权重,对地震丛集率空间分布图像分析表明,汶川M_S8.0地震前,龙门山断裂带中南段存在着长期、大范围的地震丛集率低值区,震前该段处于应力闭锁状态.对川滇地区地震丛集率低值区内背景地震与全部地震的累积次数、b值和新定义的Δb等统计参量的分析表明,龙日坝与龙门山断裂带具有地震活动的关联性,川滇地区当前的强震潜在危险区可能是巧家地区和汶川M_S8.0地震破裂尚未穿越的龙门山断裂带南段.此外,还发现b值倾向于反映局部应力场变化,而Δb能较为敏感地给出更大范围应力场的相对变化.  相似文献   

Electro‐Catalytic Degradation of Phenol Organics with SnO2‐Sb2O3/Ti Electrodes     

Kaijin Zhu  Wencai Zhang  Huifang Wang  Zhongliang Xiao 《洁净——土壤、空气、水》2008,36(1):97-102

The effect of extraordinary degradation of phenol organics on the SnO₂‐Sb₂O₃/Ti electrode is investigated through experimental research and theoretical analysis. The phenol organics contained 4‐chloro‐phenol, 4‐bromo‐phenol, and 2‐iodo‐phenol. At a current density of 4 mA cm^–2 and an electrolysis time of 12 h, the degradation efficiency of the phenols was over 98% with a relatively short degradation time, whereas the degradation time of the PbO₂/Ti electrode surpassed 40 h while delivering 100% disposal efficiency. Therefore, the effectiveness of electrochemical (EC) oxidation by the SnO₂‐Sb₂O₃/Ti was superior to that of the PbO₂/Ti electrode. At the same time, the SnO₂‐Sb₂O₃/Ti had higher oxygen generation potential and lower electron consumption than the other electrodes. This was mainly due to the effect of the middle Sb₂O₃ layer, which due to its high porosity and good catalytic effect, contributed to a better catalysis than the SnO₂ part.  相似文献   

个旧城区土壤中重金属潜在生态危害评价   总被引：28，自引：0，他引：28  

叶玉瑶  张虹鸥  谈树成 《热带地理》2004,24(1):14-17,54

采用Hakanson潜在生态危害指数法对个旧城区土壤重金属的潜在生态危害作出评价,并在此基础上应用Kriging插值方法进行了空间分析,得出了个旧城区土壤各等级生态危害的面积分布以及重金属潜在生态危害的空间变化特征.  相似文献   

二维位场在有源空间的解析延拓          下载免费PDF全文

黄宗理  郭佳 《地球物理学报》2016,59(2):704-710

位场的向下延拓对重、磁等位场资料的处理与解释具有重要意义,但现有的延拓方法在位场向下接近场源时会产生剧烈震荡,计算过程很不稳定,而且理论上都不能越过场源.本文在复数域里定义了一个具有位场和几何双重特性的新函数PFGR,提出了一套新的延拓方法,在理论上实现了位场在包含场源的下半空间里的解析延拓,经过模型试算验证,表明该方法可以延拓越过场源,准确确定场源奇点,计算过程稳定,对测量中的随机噪声也具有较好的抗干扰能力.  相似文献   

2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method          下载免费PDF全文

Yuan Yuan  Jianke Qiang  Jingtian Tang  Zhengyong Ren  Xiao Xiao 《Geophysical Prospecting》2016,64(3):767-779

To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion.  相似文献